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The exchange energy of an arbitrary collinear-spin many-body system in an exter- 
nal magnetic field is a functional of the spin-resolved charge and current densities, 
E x \n^, Tij, jjj. Within the framework of density-functional theory (DFT), we show 
that the dependence of this functional on the four densities can be fully reconstructed 
from either of two extreme limits: a fully polarized system or a completely unpolarized 
system. Reconstruction from the limit of an unpolarized system yields a generalization of 
the Oliver-Perdew spin scaling relations from spin-DFT to current-DFT. Reconstruction 
from the limit of a fully polarized system is used to derive the high-field form of the 
local-spin-density approximation to current-DFT and to magnetic-field DFT. 
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1. Introduction 

The exchange energy is a fundamental quantity of many-body physics. While it 
can be calculated analytically and exactly for simple quantum fluids, such as the 
three-dimensional interacting electron gas in the absence of magnetic fields, its 
evaluation for realistic, spatially inhomogeneous, many-electron systems is possible 
only numerically, e.g. by means of de nsity - functional theory (DFT) within the local- 
spin-density approximation (LSDA)UJ^In collinear spin-DFT (SDFT), the exact 
exchange energy is a functional of the spin-resolved charge density E x [rii, n l\- ^ n 
strong external or internal magnetic fields, orbital magnetism must be accounted 
for on t he s ame footing as spin magnetism, and SDFT is replaced by current-DFT 
(CDFTjPH or magnetic-field DFT (BDFT).EEl 

In the present paper we discuss exact and approximate expressions for the spin- 
dependence of the exchange energy E x . We show that, for arbitrary collinear-spin 
many-body systems, the dependence of E x on spin-resolved charge and current den- 
sities can be fully reconstructed from either of two extreme limits: a fully polarized 
system or a completely unpolarized system. Reconstruction from the limit of an un- 
polarized system yields a generalization of the Oliver-Perdew spin scaling relations^ 
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from SDFT to CDFT. The alternative possibility of reconstruction from the limit 
of a fully polarized system results in explict expressions for the high-field exchange 
energy, which can be used as input for the LSD A of CDFT or of BDFT. 



2. Exact properties: spin scaling 

Since exchange acts only between like-spin electrons, the exchange energy of an 
arbitrary many-electron system (homogeneous or not, in a magnetic field or not) 
can be written as the sum of two contributions, one arising entirely from the spin 
up subsystem, the other from the spin down subsystem 

E x = El + Ei (1) 

This decomposition is valid as long as there is a preferred direction relative to which 
one can define spin up and spin down, i.e. for collinear magnetism. It only ceases to 
be valid in noncollinear situations, where the quantization axis changes as a function 
of space. The exchange energy then acquires a spin off-diagonal contribution which, 
in principle, can be treated on the same footing as the diagonal termf^Here, we are 
only concerned with collinear situations, occuring in unidirectional magnetic fields. 
Under these circumstances, Eq. (fTJ always hold s 

According to the general theorems of CDFTpElthe exchange energy of an arbi- 
trary (collinear) many-body system is a functional of the spin-resolved charge and 
current densities (r), n|(r), j|(r) and jj,(r), where j(r) = Jt( r ) +j|( r ) is the para- 
magnetic current density, related to the gauge invariant physical current density 

by 

jphys(r) =j(r) + — n(r)A(r). (2) 

mc 

By symmetry, E x [n^, raj, jf, jj = E x [ni,ni, jj, jj]. The decomposition (JTJ) then im- 
plies 

E x [ni,ni,h,3i\ = E x [n h ^} +E x [n±,ii], ( 3 ) 

where the spin up and spin down contribution are given by the same functional, 
although their values for specific densities, Ej. and E^, are, generally, different. 

We now define two related functionals, the exchange energy of a spin unpolarized 
systerrP 



E» p [n,j] :=E X 



n ra j j 

2' 2' 2' 2 



(4) 



and that of a fully polarized system 

E^ p [n,j] := E x [n,0,i,0]. (5) 
By evaluating Eq. ([3]) at the densities of an unpolarized system, we obtain 

4[raJ] = i^[2ra,2j], (6) 



I 



l 
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from which it follows that we can write E x as 

E*[ni,nMl\ = \^x P [^M + \^x P [^lM> (7) 

which holds for arbitrary polarizations. By evaluating this at the densities of a 
fully polarized system, we find E x [n,j] — E FP [n,j,], which implies that the general 
exchange energy can also be written as 

£zK> n iJtJi] = E x P [ n hh] + E x P [ n i>h]- ( 8 ) 

The two limiting cases of the general functional E x [n-\, rij., jf, jjj are thus related by 

E FP [n,i} = ^ p [2n,2i}. (9) 

Equations ([7]) and ([SJ show that the exchange energy has the remarkable property 
that its functional form at arbitrary polarization is completely determined by either 
that at zero polarization or that at full polarization, without requiring interpolation 
between both limits or additional calculations inbetween them. 

If its current dependence is ignored, Eq. reduces to the spin-scaling rela- 
tions of Oliver and PerdewP which are frequently used in spin-density-functional 
theory (SDFT) to connect spin-density functionals for arbitrarily polarized systems 
to den sity funct ionals obtained from many-body calculations for an unpolarized 
systemPMJlE] Equation is the generalization of this spin scaling from SDFT 
to CDFT, and may be employed in the same way. 

Equations reconstructing the general functional from the fully polarized limit, 
such as Eq. (|8|), are less used in DFT, where one usually constructs a spin-dependent 
functional by starting from the unpolarized situation. However, they are useful in 
the special case of the high-field limit, where analytical results for the fully polarized 
exchange e nerg y of the homogeneous interacting electron gas are available for the 
first (L = O^and second (L = lj^ Landau level. 

3. Approximate expressions: high-field exchange energy of the 
homogeneous electron gas 

In the homogeneous electron gas in a constant unidirectional magnetic field B, 
the charge density and the physical current density are spatially constant. For suf- 
ficiently high magnetic fields only the lowest Landau level contributes, and the 
electrons are fully spin polarized. The per-volume exchange energy e FP (n,B) = 
E X P { n iB)./V i n this situation has been evaluated analytically by Danz and 
GlasserM^ In the limit in which all electrons occupy the bottom of the lowest Lan- 
dau level, L — 0, they find, 

e FP {n, B) = 2ire 2 X 2 n 2 [in (A 3 n) + 2.11779] , (10) 

1/2 

where A = (hc/eB) ' is the magnetic length. 

Expressions applicable to still higher fields, as well as generalization to weaker 
fields in which the restriction to the bottom of the lowest Landau level is removed 
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(while maintaining that to the lowest Lan dau level, L = 0, itself) can be found in 
the original reference In a recent papeif^l we reevaluated the excha nge en ergy 
in these limits, corrected a few minor problems in the original equations and 
discussed the range of magnetic fields and densities where each of them applies. 

For still weaker magnetic fields, two complications can take place: orbital mag- 
netism now aquires a contribution of higher Landau levels, and spin magnetism is 
reduced, as the system is not necessarily fully polarized. We dealt with the former 
complication by deriving an explicit expression for L = 1.1141 Here, we consider par- 
tial and vanishing spin polarization, which may occur for lower fields, in particular 
when the electron g factor is smaller than 2. 

In a first step, we note that the many-body calculation d^^ ^^ yield e FP (n, B) 

as function of the charge density n and the magnetic field B, i.e. the variables used 

in BDFT, but not of the spin and current densities used in CDFT. The transition 

from one set of variables to the other is affected by recalling that in the homogeneous 

electron gas in a constant magnetic field j p hys( r ) = 0, so that, from Eq. ((2]), B(n, j) = 
_r^c V x Q/ n )_ From Eq (jgj) we thus hav( p] 

Equation provides the exchange energy for arbitrarily polarized systems, as 
a function of the spin-resolved charge and current densities. As an example, we use 
it to generate from Eq. (p~0|) the exchange energy of an unpolarized system. This can 
be done either by using pT|) to generate the result for an arbitrarily polarized system 
and subsequently substitute the densities of an unpolarized system, or, equivalcntly, 
directly from the special case of Eq. © . The result is 



e 



NP 



(n, B) = 7re 2 A 2 n 2 [in (\ 3 n) + 1.42464] . (12) 

All other relations, valid for different magnetic field ranges can be treated 
in the same way. We have chosen to exemplify the procedure for Eq. ([TO]) . because 
this allows an indepen den t numerical test of our previous assertion (which was based 
on explicit calculation^!) that the Danz-Glasser expression (flTJ|) is correct, and not 
the alternative resultP^l 

e F p ' alt (n, B) = 27re 2 A 2 n 2 [in (\ 3 n) + 2.32918] , (13) 

which, when scaled to the unpolarized case, becomes 

e% p ' alt (n, B) = 7re 2 A 2 n 2 [in (\ 3 n) + 1.63603] . (14) 

I -j hyl 

In the unpolarized case, Takada and Goto numerically calculated e x , so that we 
can use their data to check the validity of Eqs. (Tj"2")) and (Ti"4"|) , and thus also of (JT0J) 
and (Ti"3"|) . We find that (fT2")) is in perfect agreement with numerical data extracted 
from figure 2(b) of Takada and Goto whereas (fT4"]) deviates from the numerical 
data by ~ 10% in the metallic density range. 
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4. Conclusions 

Two types of spin-scaling relations have been demonstrated for the current- 
dependent exchange energy, one employing the fully polarized limit and the other 
the unpolarized limit. Simple connections, valid for arbitrary many-body systems, 
exist between these limits and the general functional E x [n^ , n±, jj, jj. In the spe- 
cial case of the homogeneous electron g as , such c onnections , applied to many-body 
calculations in the fully polarized limit JHHHHI yield explicit formulas that allow 
to use available numerical data for unpolarized systems^ to arbitrate between con- 
flicting expressions at full polarization. Results for arbitrary polarization can be 
used as input for the local-spin-density approximation to current-density-functional 
theory. 
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